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Abstract 

This  thesis  presents  a  career  model  of  the  USAF  Pave  Low 
helicopter  pilot  force.  A  SLAM  II  computer  simulation  was 
developed  to  permit  experimentation  with  key  decision 
variables.  Response  surface  methodology  was  used  to  produce 
two  response  equations  which,  In  turn,  became  goals  for  a 
goal  programming  model.  This  model  provided  optimal  values 
for  the  decision  variables  In  order  to  achieve  manning  and 
experience- level  goals. 


MODELING  THE  PAVE  LOW 
HELICOPTER  PILOT  CAREER  FIELD 


I .  Introduct.ion 


Tha  Issue 

“Low-lntenslty  conflict  represents  the  most  probable 
arena  for  the  application  of  O.S.  military  force  in  the 
foreseeable  future,"  according  to  Major  Kenneth  Page  in  his 
article,  "OS  Air  Force  Special  Operations:  Charting  a 
Course  for  the  Future."  As  a  result.  Congress  has  directed 
the  OSAF  to  improve  its  special  operations  helicopter 
capability  by  building  up  its  small  force  of  MH-53  Pave  Low 
helicopters.  (Page,  1987:58,66) 

However,  at  the  same  time,  many  other  OSAF  helicopters 
are  being  phased  out  of  service,  resources  which  in  the  past 
have  provided  flying  hours  and  experience  for  pilots  who 
would  later  transition  to  the  more  sophisticated  Pave  Low. 
Therefore,  there  is  concern  that  the  Air  Force  may 
eventually  be  faced  with  the  prospect  of  sending  primarily 
new,  less  experienced  pilots  to  the  Pave  Low,  pilots  who 
have  not  had  the  benefit  of  developing  airmanship  skills  in 
a  less  complex  helicopter. 

Problem  Statement 

As  the  Air  Force's  pool  of  experienced  helicopter  pilots 
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shrinks,  sustaining  an  experienced  Pave  Low  pilot  force  will 
become  more  difficult,  AFHPC  does  not  have  a  model  of  the 
Pave  Low  pilot  force  to  study  force  sustainment,  so  they  are 
not  able  to  predict  how  policies  will  affect  the  force. 
Baaaarch..  Qbj active  and  Scope 

The  objective  of  this  research  is  to  develop  and  analyze 
a  model  of  the  1025N  (Pave  Low  pilot)  career  field  as  it 
will  exist  after  the  force  reaches  full  strength  in  order  to 
examine  how  to  best  sustain  the  pilot  force  in  the  future. 

This  investigation  will  be  limited  to  Pave  Low  pilots  and 
copilots,  from  the  time  they  enter  the  Pave  Low  school  until 
they  stop  flying  the  Pave  Low  permanently.  It  will  include 
the  time  those  individuals  spend  in  career-broadening/staff 
assignments,  assuming  they  return  to  the  Pave  Low  cockpit. 

It  will  not  include  other  helicopter  pilots  except  as  inputs 
to  the  Pave  Low  force. 

Key  Terms /Acronyms 

AFMPC  -  Air  Force  Military  Personnel  Center 

MAC  -  Military  Airlift  Command 

PME  -  Professional  Military  Education 

POP  -  Prior  Qualified  Pilot 

Rated  Officer  -  a  pilot  or  navigator 

Rated  Supplement  -  rated  officers  holding  non-flying  jobs  in 
support  career  fields 

OHT  -  Undergraduate  Helicopter  Training  or  a  recent  graduate 
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Chapter  II  provides  a  background  to  rated  officer 
management  and  reviews  modeling  approaches  used  in  the 
past  to  model  personnel  systems.  Chapter  III  covers  the 
modeling  approach  and  output  analysis  used  in  this  research. 
Chapter  IV  is  the  results  of  the  analysis  and  Chapter  V 
contains  the  conclusions  drawn  from  that  analysis. 
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II. 


Introduction 

This  chapter  describes  how  the  Air  Force  manages  rated 
officers  and  then  gives  brief  descriptions  of  various 
methodologies  that  have  been  used  to  model  personnel 
systems:  In  particular,  military  personnel  systems.  This 

review  will  help  In  the  selection  of  the  most  appropriate 
modeling  technique  for  the  Pave  Low  system.  The  literature 
Is  reviewed  In  a  topical  order:  rated  officer  management 
and  personnel  modeling  approaches,  which  Include  both 
prescriptive  and  descriptive  models. 

Bated  Officer  Management 

AFR  36-20  describes  the  gate  system  the  Air  Force  uses 
to  manage  Its  rated  officer  force.  Rated  officers  are 
required  to  complete  six  years  of  operational  flying  prior 
to  being  assigned  to  non-flying  duties.  Furthermore,  "it  is 
Air  Force  policy  that  as  many  members  as  possible  perform  at 
least  9  years  of  operational  flying  duty  during  the  first  18 
years  of  aviation  service"  (AFR  36-20:11). 

Rated  resource  managers  at  AFMPC  who  manage  the 
assignments  of  rated  officers  roust  balance  these  Air  Force 
requirements  with  the  need  to  provide  officers  the 
opportunity  to  broaden  their  careers  in  non- flying 
assignments  such  as  staff  duty,  rated  supplement,  AFIT, 
and  PME. 
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Additionally,  the  resource  manager  must  consider  what 
impact  moving  experienced  personnel  has  on  the  flying 
organisation.  If  a  less-experienced  replacement  pilot  must 
be  absorbed  into  the  unit,  he  must  be  allowed  to  fly  to  gain 
sufficient  experience  so  he  becomes  an  asset  to  the 
organisation.  Absorption  capacity  is  "currently  the  most 
Important  factor"  in  determining  the  production  rates  of 
undergraduate  flying  training  (Rated  Management  Document, 
1987:6-1) . 

Davie  (1986)  said  absorption  could  be  a  major  problem  for 
the  Pave  Low  force  as  it  expands  and  the  pool  of  experienced 
helicopter  pilots  shrinks.  Less-experienced  pilots  must 
have  time  to  develop  basic  flying  skills  in  addition  to 
learning  the  Pave  Low  aircraft  atnd  mission.  A  linear 
programming  model  can  be  used  to  demonstrate  how  experience 
and  stability  requirements  can  severely  restrict  the  number 
of  inexperienced  pilots  (OHTs)  the  Pave  Low  force  can  absorb 
and  still  maintain  its  experience  standard  (Rated  Management 
Document,  1987:6-3). 

AFMPC  recognizes  that  while  a  linear  programming  model 
can  Indicate  an  optimal  mix  of  experienced  and  inexperienced 
pilots  as  Inputs  to  the  Pave  Low  force,  that  optimal  mix 
'Tflll  not  necessarily  reflect  what  is  available  for 
assignment  (Rated  Management  Document,  1987:6-5).  AFMPC 
needs  a  model  of  the  force  that  will  allow  it  to  vary,  among 
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other  things,  the  mix  of  pilots  to  study  the  effects  on 
force  sustainment. 

Personnel  Modeling  Approaches 

Personnel  models  in  the  literature  generally  fall  into 
one  of  two  broad  categories:  prescriptive  or  descriptive. 
Prescriptive  models  prescribe  an  optimal  solution  fo^  a 
given  set  of  conditions.  They  generally  involve  some  type 
of  math  programming  such  as  linear  programming  (LP),  goal 
programming,  or  network  flow  programming.  Descriptive 
models  are  meant  to  imitate  the  systems  and  provide  insight 
about  that  system.  A  common  descriptive  model  is  the 
simulation  model  (Olson,  1987:14,17). 

Prescriptive  Models. 

Linear  Programming.  A  linear  programming  model  is  a 
set  of  linear  equations  that  bounds  a  feasible  region.  An 
objective  function  is  used  to  determine  a  point  within  that 
region  that  represents  an  optimal  solution.  The  simplex 
method,  developed  by  Dantzlg  (Hillier  and  Lleberman, 

1986:4),  can  be  used  to  solve  a  linear  programming  problem. 

Klelnman  and  Qoudreau  (1977)  used  linear  programming  to 
determine  the  optimal  number  of  Havy  officer  accessions  from 
nine  commissioning  sources.  Within  the  constraint  to  fill 
required  slots,  the  model  minimizes  life  cycle  costs. 

Aronson  and  Thomson  (1985)  applied  forward  simplex  to  the 
multi-period  linear  goal  programming  problems  of  Charnes, 
Cooper,  and  Niehaus  (Charnes  et  al.,  1972).  They  used 
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Markov  transition  probabilities  between  personnel  grades  to 
determine  prooiotlons  from  one  period  to  the  next. 

TOPOPS  (Akam  and  Nordhauser,  1974)  Is  a  OSAF  officer 
procurement  program  that  determines  the  best  combination  of 
officer  recruits  from  the  various  commissioning  sources.  It 
can  be  used  to  minimize  costs  or  to  maximize  the  quality  of 
20  officer  types  over  a  five  year  period.  Charpie  (1987) 
used  both  linear  programming  and  simulation  to  study  the 
DSAF  B-52  navigator  career  field.  The  linear  program 
determined  the  optimal  state  of  the  force,  Including  optimal 
accession  and  training  rates  for  a  given  scenario.  As  a 
basis  for  his  LP  model,  Charpie  used  a  cross-sectional 
network  model  such  as  the  one  In  figure  1 . 


State  Tisa  t  Tiae  t^-l 
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Figure  1 .  Cross-Sectional  Network  Flows 
(Charpie,  1987) 
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This  network  represents  positions  or  job  types  as  nodes  with 
I  ^ 

the  connecting  arcs  representing  possible  movement  from  one 
job  (at  time  t)  to  another  (at  time  t+1) . 

Charpie's  LP  model  has  five  categories  of  constraints* 

I 

1)  conservation  of  flow,  2)  manning  level  requirements  for 
the  various  jobs,  3)  minimum  and  maximum  flows  from  one 
position  to  another,  4)  attrition,  and  5)  tour  length.  The 

I 

objective  function  seeks  to  minimize  a  variety  of  variables 
including  overall  manning  requirements,  attrition, 
accessions,  and  upgrade  training. 

I 

Charple  used  the  LP  model  to  establish  an  optimal  steady 
state  condition  for  the  navigator  career  field  for  a  given 
scenario.  The  model,  however,  did  have  drawbacks.  Charpie 
reported  that  the  LP  model  could  not  incorporate  as  many 
different  aspects  of  the  career  field  as  the  simulation 
model  could,  it  tended  to  Ignore  unrealistic  personnel 
allocations,  and  finally,  the  solution  was  merely  a 
"snapshot"  in  time  and  could  not  reflect,  for  example, 
attrition  over  a  period  of  time. 

Goal  Programming-  Goal  Programming  is  an  adaptation 
of  linear  programming  that  includes  multiple  goals  in  the 
objective  function.  Chames,  Cooper,  and  Niehaus  (1968) 
developed  several  different  personnel  models  for  the  DS  Navy 
using  goal  programming  embedded  with  Markov  processes. 

Siverd  and  Thomson  (1979)  used  ratio  goal  programming  to 
assign  personnel  with  various  skill  levels  to  jobs  based  on 
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a  preferred  olx  of  those  skill  levels.  By  doing  so,  they 

I 

were  able  to  maximize  a  measure  of  organizational 
effectiveness . 

Network  Flow  Prngramming.  Network  models  can  be 

I 

visualized  as  a  set  of  nodes  connected  by  arcs  or  paths.  In 
a  typical  personnel  network  model,  the  nodes  could  represent 
jobs  with  individual  workers  flowing  along  the  arcs,  from 

I 

node  to  node,  as  they  change  jobs.  Thompson  (1976)  used  a 
network  transshipment  model  to  approximate  a  linear 
programming  model  of  a  personnel  system. 

I 

O'Conner  (1982)  developed  an  interactive  computer  program 
that  used  networks  to  model  the  Navy  pilot  force.  The  model 
was  used  to  test  the  Implications  of  various  policy 

I 

alternatives  on  Navy  pilot  requirements. 

Liang  and  Thompson  (1986)  developed  an  assignment 
optimization  model  that  uses  a  capacitated  network  to  match 
Navy  enlisted  personnel  to  jobs  according  to  Navy  policy. 
People  and  jobs  are  both  modeled  as  nodes  with  a  connecting 
arc  if  an  individual  is  eligible  for  that  particular  job. 

Olson  (1987)  used  a  single-commodity  network  flow  with 
side  constraints  to  model  the  DSAF  rated  officer  force.  The 
model  balanced  flying  and  non-flying  assignments  to  ensure 
most  rated  officers  met  their  flying  gates. 

Olson's  network  is  four-dimensional  in  the  following 
parameters:  1)  time  period  under  examination,  2)  duty 

assignment,  3)  aviation  service  date  (ASD)  year  group,  and 
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4)  accumulated  flying  gate  credit.  Each  node  represents 
time  (In  years),  the  first  parameter,  and  an  Individual's 
status ,  made  up  of  the  other  three  parameters .  The  arcs 
between  the  nodes  represent  assignments  an  officer  could 
serve  to  get  from  one  node/status  to  another  node/status . 
Olson's  network  contains  975  nodes  and  2374  allowable  arcs. 

The  three  dimensional  network  In  figure  2  was  taken 
directly  from  Olson's  thesis.  It  Illustrates  how  three  of 
the  four  parameters  Interact.  Arc  1  represents  a  one  year 
assignment  to  duty  type  A,  with  a  resulting  Increase  In 
flying  gate  credit.  Arc  2  represents  a  two  year  assignment 
to  duty  type  B  with  no  Increase  In  flying  gate  credit. 

The  parameter  not  depicted  In  figure  2,  ASD  year  group. 
Increases  directly  with  time.  Thus  the  ASD  year  group  for 
the  Individual  at  node  B3/I  Is  two  years  greater  than  It  was 
when  he  was  at  node  Al/I . 

In  addition  to  network  flow  constraints  that  prevent 
Illogical  assignments  (such  as  AFIT  followed  by  AFIT),  Olson 
used  side  constraints  to  ensure  that  the  model  assigns 
sufficient  personnel  to  each  duty  type  without  overmanning 
any  of  them.  Side  constraints  were  also  used  to  ensure 
units  are  manned  by  sufficiently  experienced  personnel . 

Olson  used  an  optimization  routine  called  NETSID  to  solve 
his  formulation.  NETSID  Is  a  simplex  linear  programming 
algorithm  that  can  efficiently  Incorporate  side  constraints. 
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Figure  2.  Three-Dimensional  Assignment  Network 

(Olson,  1987) 


Olson  defines  a  side  constraint  as  any  constraint  that 
applies  to  flow  across  multiple  arcs.  The  problem  was 
formulated  as: 

Minimize  cx 
Subject  to:  Ax  =  r 


Sx  =  b 
i. 

0  i  X  i.  u 


11 


A  Is  a  matrix  with  a  1  if  the  arc  flows  out  of  the  node, 
-1  if  the  arc  flows  into  the  node,  and  0  otherwise.  S  is  a 
matrix  of  side  constraints  with  b  as  the  right  hand  side. 

An  example  of  a  side  constraint  would  be  one  that  required  a 
minimum  of  60  percent  of  all  flying  duty  positions  to  be 
occupied  by  pilots  in  ASD  year  group  six  or  higher.  To 
determine  which  arcs  are  affected,  the  program  must  consider 
all  arcs  terminating  at  flying  duty  nodes  with  an  ASD  year 
group  greater  than  or  equal  to  six.  The  total  of  these 
flows  would  have  to  be  greater  than  or  equal  to  60  percent 
of  the  total  flow  into  all  flying  duty  nodes. 

The  r  vector  represents  the  requirements  at  each  node  and 
u  is  the  vector  of  upper  bounds  for  the  arc  flows .  The  x 
vector  is  the  solution  vector.  The  objective  function 
minimizes  total  costs  associated  with  failing  to  allow 
pilots  to  meet  their  flying  gates  by  the  time  they  get  to  a 
given  ASD  year  group.  The  cost  vector,  c,  represents  the 
costs  along  arcs  that  represents  such  failures .  For 
example,  any  node  with  an  ASD  value  of  12  years  or  more  and 
a  gate  credit  value  of  fewer  than  six  years  represents  a 
failure  to  meet  the  first  (six  year)  gate.  Any  arc 
terminating  at  such  a  node  would  have  a  penalty  or 
cost  associated  with  it. 

Paacriptive  Models- 

Simulation  Models.  Computer  simulations  of  personnel 
systems  generally  model  at  either  the  individual  (entity)  or 
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the  group  (aggregate)  level.  Entity  models  provide  more 
I  detail,  but  can  also  require  more  input  data  and  longer 

computer  run  times  (Clark  and  Lawson,  1984:42).  Leupp 
(1969)  describes  an  entity  simulation  model  the  U.S.  Havy 
I  uses  to  project  the  number  of  SEABEE  personnel  needed  in  the 

future.  The  model  considers  accessions,  attrition, 
promotions,  and  assignment  rotations. 

\  Looper  (1979)  describes  a  large-scale  entity  simulation 

model  called  the  Career  Area  Rotational  Model  (CAROM)  that 
uses  Monte  Carlo  techniques  to  model  the  career  progression 
)  of  airmen  in  the  OSAF.  The  model  updates  accessions, 

promotions,  assignment  rotations,  and  attrition  on  a  monthly 
basis  within  a  single  career  field  for  up  to  30  years. 

'  Percich  (1987)  simulated  the  career  progression  of  OSAF 

strategic  airlift  pilots  as  a  basis  to  study  moving  costs. 
The  model  is  a  network  and  discrete-event  simulation.  He 

*  built  the  network  representation  of  the  career  field  using 
the  SLAM  II  language  (Pritsker,  1986).  The  entities  that 
flow  through  the  system  are  the  individual  pilots,  modeled 

*  from  the  time  they  enter  undergraduate  pilot  training  (UPT) 
until  they  leave  the  Air  Force.  As  each  pilot  moves  from 
location  to  location  with  reassignments,  the  model  collects 
data  on  moving  costs  associated  with  each  pilot. 

While  the  network  structure  of  SLAM  II  is  adequate  to 
model  the  flow  of  pilots  from  one  assignment  to  another, 
Percich  wrote  separate  subroutines  to  initialize  the 
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network,  set  tour  lengths,  model  attrition,  make 
assignments,  and  collect  moving  costs. 

Using  the  results  of  experimentation  with  the  simulation 
model,  Perclch  developed  two  response  surfaces.  He  first 
group  screened  a  total  of  12  predictor  variables  or  factors, 
divided  Into  five  homogeneous  groups.  Using  linear 
regression,  he  was  able  to  Identify  three  significant  groups 
containing  a  total  of  eight  factors.  He  then  screened  those 
remaining  eight  factors  down  to  four  significant  factors. 
Based  on  those  four  factors,  he  built  both  linear  and 
quadratic  metamodels,  polynomials  that  are  significantly 
easier  and  quicker  to  use  than  a  simulation  model . 

Charpie  (1987)  developed  a  model  of  the  B-52  navigator 
career  field  In  the  U.S.  Air  Force  and  used  It  to  examine 
changes  In  policies  affecting  the  crew  force.  He  used  a 
network  representation  In  which  nodes  represent  the  various 
Jobs  that  a  navigator  can  hold  with  navigators  as  the 
entitles.  The  arcs  represent  allowable  assignments  for 
Individuals  at  each  stage  of  their  careers. 

Lltko  and  Travis  (1982)  developed  a  dynamic  simulation 
model  to  examine  the  potential  consolidation  of  four  USAF 
career  fields  Into  one.  The  model  used  a  combination  of  the 
network,  discrete-event,  and  continuous  flow  capabilities  of 
SLAM.  The  flow  of  personnel  to  and  from  overseas 
assignments  was  modeled  as  a  continuous  process,  while  the 
actual  consolidation  of  the  career  field  was  modeled  as  a 
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network  activity.  The  discrete-event  portion  of  the  model 
was  used  to  subtract  or  add  personnel  to  the  career  field  as 
technical  school  classes  enter  and  graduate.  The  school 
takes  personnel  out  of  the  system  for  the  duration  of  the 
school . 

Perclch  (1987)  reports  the  Integrated  Simulation 
Evaluation  Model  (Prototype)  (ISEM-P)  is  a  large-scale 
aggregate  simulation  that  AFMPC  can  use  to  model  the  entire 
OSAF  personnel  system.  Starting  with  a  mission  plan,  it 
establishes  manpower  requirements  and,  after  considering 
attrition,  determines  the  necessary  accessions  to  build  the 
required  inventory. 

Forrester  (1980)  developed  systems  dynamics  to  incorporate 
closed  feedback  loops  into  simulation  models  in  an  effort  to 
more  realistically  represent  real  world  interactions. 

Knight  (1978)  used  systems  dynamics  in  a  highly  simplified, 
aggregate  model  of  OSAF  pilots.  He  used  the  model  to 
examine  policies  controlling  the  size  of  undergraduate  pilot 
training  classes  and  Instructor  to  student  ratios. 

Lawson  (1982)  used  an  aggregate  simulation  model  of 
enlisted  force  structures  within  A.F.  Communications  Command 
(AFCC)  to  experiment  with  policies  regarding  the  flow  and 
distribution  of  personnel .  His  model  also  employed  the 
system  dynamics  approach. 

Clark  and  Lawson  (1984)  used  the  systems  dynamics 
approach  to  model  a  segment  of  the  USAF  enlisted  force  that 
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is  often  assigned  overseas .  They  started  with  an  influence 
diagram  that  graphically  showed  both  feedback  structures  and 
cause  and  effect  relationships.  Their  simulation  model 
required  significantly  less  computer  code  and  structure  than 
it  would  have  if  they  had  used  an  entity  model. 

The  models  described  in  this  chapter  provide  a  basis  on 
which  to  choose  a  methodology  to  model  the  Pave  Low  pilot 
force.  The  next  chapter  describes  that  methodology. 
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III. 


The  purpose  of  this  chapter  is  to  examine  the  data 
available  to  build  a  model  and  then  to  discuss  the  modeling 
approach  and  the  model  Itself.  Included  are  discussions  of 
the  model  assumptions,  verification,  validation,  and 
experimental  design. 


Because  the  steady  state  Pave  Low  personnel  system  to  be 
modeled  Is  not  yet  a  reality,  data  on  the  system  Is  scarce. 
The  relatively  new  career  field  has  tended  to  hold  pilots  in 
the  cockpit  longer  than  in  other  weapon  systems  due  to  the 
complex  nature  of  the  mission  and  equipment.  Furthermore, 
there  is  currently  only  one  assignment  location,  whereas 
there  will  eventually  be  multiple  locations.  Because  this 
research  is  intended  to  examine  a  steady  state  environment, 
data  on  average  assignment/ tour  lengths,  number  of  tours, 
and  flying  hours  taken  from  the  current  system  would  be 
misleading  at  best. 

As  a  result,  most  of  the  numbers  used  in  this  research 
have  been  gleaned  from  discussions  with  individuals  at 
AFHPC,  MAC,  and  23rd  AF.  They  are  educated  guesses  at  what 
the  steady  state  system  will  look  like.  The  exception  is 
the  retention  data  (see  appendix  A),  which  is  actual 
historical  data  reflecting  the  entire  helicopter  force  for 
the  year  1987. 
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For  the  purpose  of  this  investigation,  it  will  be  assumed 
that  the  system  (career  field)  to  be  modeled  is  in  steady 
state.  That  is,  all  the  Pave  Low  helicopters  are 
operational  and  in  place,  the  units  are  fully  manned,  and 
some  percentage  of  the  pilots  are  temporarily  in  non-flying 
assignments.  Since  this  steady-state  system  does  not  yet 
exist,  however,  the  model  will  use  reasonable  estimates  for 
parameters  such  as  average  time  to  upgrade  from  one  crew 
position  to  another. 

A  sustained  pilot  force  will  bo  defined  as  100  fully 
qualified  pilots  (not  including  those  in  staff  positions), 
60%  of  whom  must  be  experienced  Pave  Low  pilots.  An 
experienced  Pave  Low  pilot  is  one  with  at  least  850  total 
flying  hours,  300  of  which  must  be  in  the  Pave  Low 
helicopter  (AFR  51-2).  Finally,  these  100  pilots  are 
assumed  to  be  distributed  across  the  different  crew 
positions  and  ranks  as  they  would  in  a  real  world  flying 
squadron . 

As  previously  mentioned,  retention  data  used  is  from 
historical  data.  This  data  remains  constant  throughout  the 
simulation,  although  in  reality  these  figures  would  change 
somewhat  from  year  to  year.  In  fact,  all  system  parameters 
such  as  input  rates  into  the  system  (INP),  the  mix  of  UHTs 
(OHT)  versus  PQPs,  and  tour  lengths  (TRL)  are  constant 
throughout  the  simulation.  In  the  real  world,  these 


necessarily  fluctuate  with  changlns  requirements  and 
policies . 

Promotions  to  higher  grades  and  upgrades  to  higher  crew 
positions  are  not  modeled.  It  Is  assumed  these  occur  at 
normal  rates  and  times.  In  reality.  If  It  were  perceived 
that  promotions  and  upgrades  lagged  behind  those  In  other 
career  fields,  retention  rates  would  probably  decline. 
Modeling  Approach  Selection 

The  models  discussed  so  far  can  be  categorized  as  dynamic 
or  static,  stochastic  or  deterministic,  aggregate  or  entity, 
and  proscriptive  (optimization)  or  descriptive  (non- 
optlmlzatlon) .  The  relative  advantages  and  disadvantages  of 
these  various  approaches  depend  on  the  system  to  be  modeled. 

Personnel  systems  are  by  nature  dynamic,  that  is, 
changing  over  time.  Static  models  generally  ignore  the 
passage  of  time.  Instead  taking  a  "snapshot”  In  time  of  that 
system.  Lltko  (1982:33)  says  this  method  can  be  misleading 
when  modeling  a  changing  process  such  as  a  personnel  career 
field.  For  this  reason,  a  dynamic  model  seems  to  be 
appropriate  for  this  research. 

Similarly,  real-world  personnel  systems  have  a  degree  of 
unpredictability  about  them.  To  treat  a  career  field  as 
deterministic  suggests  the  ability  to  predict  or 
predetermine  factors  such  as  tour  lengths  or  individual 
career  decisions.  A  stochastic  model,  on  the  other  hand, 
uses  random  numbers  to  model  the  variability  In  such  factors 


19 


and  will  therefore  be  used. 

Lltko  (1982:36)  says  the  disadvantages  of  an  entity  model 
for  large  personnel  systems  are  the  need  for  an  extensive 
data  base  and  longer  computer  run  times.  The  disadvantage 
of  the  aggregate  approach  Is  the  potential  loss  of  necessary 
detail  and  validity. 

The  Pave  Low  pilot  data  base  Is  currently  very  small 
since  It  Is  a  relatively  new  weapon  system.  Even  when  It 
reaches  Its  projected  peak,  It  will  consist  of  just  over  100 
actively  flying  pilots  with  about  that  many  In  non-flying 
positions .  With  so  few  entitles  to  model ,  the  need  for 
detail  seems  to  outweigh  the  difficulty  of  dealing  with  a 
large  data  base.  Furthermore,  Initial  runs  of  a  simple 
prototype  entity  simulation  model  suggest  very  short 
computer  run  times:  only  a  few  seconds  actual  CPU  time. 

The  question  of  whether  to  use  an  optimization  or  non- 
optlmlzatlon  model  should  be  tied  to  the  type  of  questions 
to  be  answered  with  the  model .  The  primary  objective  of 
this  research  is  to  provide  AFMPC  with  a  tool  to  study  the 
impact  (on  the  Pave  Low  pilot  force)  of  policy  decisions, 
particularly  decisions  concerning  the  allowable  mix  of  UHTs 
versus  PQPs  as  Inputs  to  the  Pave  Low  system.  This  research 
is  Intended  to  provide  "an  array  of  acceptable  policies  . . . 
from  which  the  planner  can  choose  the  one  with  the  best 
tradeoff."  (Jaquette  et  al . ,  1977:6).  Furthermore,  it 
should  allow  decision  makers  to  ask  "what  If"  questions  and 
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evaluate  tentative  results  without  actually  Implementing  any 
policies.  A  secondary  objective  Is  to  determine,  for  a 
given  set  of  circumstances,  an  optimal  policy. 


To  meet  the  first  objective,  the  modeling  approach  of 
choice  will  be  a  dsmamlc,  stochastic,  entity  simulation 
(non-optlmlzatlon)  model.  Response  surface  methodology  will 
be  applied  to  the  simulation  output  to  develop  response 
equations  that  describe  the  responses  In  terms  of  policy 
variables.  These  equations  will.  In  turn,  be  used  as  goals 
In  a  goal  programming  (optimization)  model  to  determine  an 
optimal  policy,  one  that  meets  the  manning  and  experience 
goals  as  closely  as  possible. 
tlQdftl.  OyarvieH 

The  Pave  Low  simulation  model  was  fashioned  after 
Percich's  Strategic  Airlift  PCS  cost  simulation  model 
(Perclch,  1987).  Perclch  modeled  a  rated  career  field,  the 
strategic  airlift  pilots.  In  order  to  study  permanent  change 
of  station  (PCS)  costs.  His  model  Is  an  entity  simulation 
model  written  In  SLAM  II,  a  language  which  Is  also  available 
for  this  research.  Like  that  model,  the  Pave  Low  model 
breaks  the  pilot's  career  Into  three  periods:  one  to  six, 
seven  to  eleven,  and  twelve  to  seventeen  years.  The  period 
from  18  years  to  retirement  Is  modeled,  but  not  considered 
for  flying  assignments. 

The  first  period,  from  one  to  six  years,  consists  of 
flight  school  for  OHTs,  Pave  Low  school  (CCTS),  and  the 
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first  fl/lng  assignmentCs) .  Hon-f lying  assignments  are  not 
allowed  during  this  period  In  keeping  with  gate 
requirements.  Subsequent  periods  Include  both  flying  and 
non-flying  assignments.  Figure  3  Is  a  macro  flow  chart  of 
the  Pave  Low  Simulation  Model . 

Intended  Dse  of  the  Model 

The  Pave  Low  simulation  model  Is  Intended  to  provide 
decision  makers  a  range  of  options  to  sustain  the  Pave  Low 
force  at  required  manning  and  experience  levels.  It  allows 
decision  makers  to  ask  “what  If”  questions;  that  Is,  to 
change  key  parameters  and  see  the  potential  Impact  a  policy 
decision  might  have  on  the  actual  system. 

Pave  Low  Simulation  Model 

The  model  Itself  Is  written  In  SLAM  II.  SLAM  II  Is  a 
FORTRAN  based  modeling  language  with  the  capability  to 
Incorporate  both  discrete-event  and  continuous  modeling  Into 
a  network  model  (Prltsker,  1984). 

SLAM  II  Network  Code.  Entity  attributes,  values 
associated  with  particular  entitles,  and  global  variables, 
which  are  not  associated  with  any  particular  entity,  are 
defined  In  table  1 .  A  network  diagram  and  the  network 
computer  code  appear  In  appendices  B  and  C.  The  first 
create  node  creates  an  entity  whose  sole  purpose  Is  to  call 
subroutine  CLEAR  at  200  time  units  (months).  Subroutine 
CLEAR  clears  all  statistical  arrays  of  transition  period 
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Figure  3.  Pave  Low  Model  Flow  Chart 
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data  in  order  to  eliminate  bias  created  by  the  warm-up 
period  prior  to  steady  state . 


The  second  create  node  creates  an  entity  every  XX(6) 
months  that  represents  a  CCTS  class  of  two  student  pilots. 


Table  1  SLAM  II  Attributes  and  Variables 


Attribute 


1 

tour  length 

2 

hours  per  tour 

3 

pilot  type  (OHT  or  PQP) 

4 

time  in  service 

5 

Pave  Low  flying  hours 

6 

hours  until  experienced 

7 

separation  indicator 

8 

Global 

Variable 

1 

number  of  experienced  pilots 

2 

number  of  actively  flying  pilots  =  TOT 

3 

hours  required 

4 

XX(1)/XX(2)  =  EXP 

5 

excess  hours 

6 

time  between  class  creations 

7 

percent  OHT 

8 

low  side  tour  length  distribution 

9 

mode  of  tour  length  distribution 

10 

high  side  tour  length  distribution 

11 

hours  payback  (PAY) 

_ 12- 

XX(7)  percent  of  those  pilots  are  labeled  UHTs  and  have  a 
time-in-service  of  12  months  (for  flight  school  and  a  PCS). 
The  rest  of  the  pilots  are  labeled  PQPs  and  have  a  time  in 
service  drawn  from  a  triangular  distribution  with  low,  mode, 
euxd  high  values  of  42,  48,  and  54  months  respectively.  In 
the  absense  of  a  known  distribution,  the  trlagular 
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distribution  was  used  to  provide  a  range  of  possible  values 
with  the  mode  being  the  most  likely  value.  This  represents 
approximately  a  year  for  flight  school  and  a  three  year  tour 
In  another  type  of  helicopter.  All  students  are  sent  to 
CCT£  for  seven  months  (accumulating  80  Pave  Low  hours),  at 
which  time  they  are  counted  as  part  of  the  actively  flying 
line  force. 

Hext,  the  pilots  are  assigned  a  tour  length  drawn  from  a 
triangular  distribution,  again,  since  It  Is  not  a  known 
distribution.  The  event  node  calls  subroutine  FLIER,  which 
determines  retention/separation  status  and  updates 
applicable  attributes,  such  as  Pave  Low  flying  hours.  If  a 
DHT  needs  more  flying  time  In  order  to  reach  850  hours,  he 
Is  looped  back  for  a  second  assignment. 

The  next  section  of  the  network  determines,  to  the  month, 
when  each  pilot  becomes  experienced  and  immediately  updates 
the  ratio  of  experienced  to  total  active  pilots.  Pilots  who 
separate  before  becoming  experienced  are  collected  at  a 
different  node  from  those  who  separate  after  becoming 
experienced . 

Pilots  who  have  not  finished  paying  back  for  their  Pave 
Low  training  are  looped  back  for  another  flying  assignment. 
The  rest  move  on  to  the  7  to  11,  or  the  12  to  17  year  group, 
as  appropriate.  From  Queue3,  XX(12)/2  pilots  are  blocked 
for  flying  assignments  beginning  sometime  during  the  7  to  11 
year  period.  The  same  number,  XX(12)/2  are  blocked  in  the 
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12  to  18  year  group.  XX(12)  is  the  total  number  of  blocked 
cockpits,  but  these  are  divided  between  two  career  stages. 
After  these  slots  are  filled,  the  remainder  of  the  pilots 
fill  non-flying  assignments.  These  non-flying  pilots  are 
subtracted  from  the  force  and  the  experience  level  Is 
recalculated. 

The  12  to  17  year  group  part  of  the  network  parallels  the 
previous  section.  The  18  year  to  retirement  section  was 
only  used  for  non-flying  assignments  for  this  research,  but 
could  be  used  to  block  cockpits  for  more  senior  pilots.  The 
final  section  of  the  model  Is  data  collection  nodes. 

gQRTBAW  Subroutines.  FORTRAN  subroutines  are  located  in 
appendix  D.  Subroutine  INTLC  Is  used  to  read  In  the 
parameter  values  from  a  coded  design  matrix.  It  then 
decodes  the  parameters  Into  a  form  meaningful  to  the  SLAM  II 
network.  The  subroutine  also  holds  the  retention  data 
table.  Subroutine  INTLC  could  be  used  to  Initialize  or  load 
the  network  to  shorten  the  transition  period  to  steady 
state . 

Subroutine  FLIER  rounds  the  time  ’In  service  and  pending 
tour  length  to  whole  years  in  order  to  apply  the  retention 
data.  It  then  makes  a  draw  from  a  uniform  (0,1) 
distribution  and  compares  that  number  to  the  applicable 
retention  figures.  If  the  random  draw  Is  greater  than  the 
retention  figure  for  any  year  of  the  pending  assignment,  the 
tour  length  Is  shortened  and  attributes  are  updated 
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accordingly.  The  entity  Is  given  a  separation  indicator 
that  will  cause  it  to  leave  the  network.  If  the  entity  does 
not  separate,  flying  hours  2uid  time-ln-servlce  are  updated 
according  to  the  originally  assigned  tour  length. 

Subroutine  NOHFLY  is  Identical  to  subroutine  FLIER,  except 
flying  hours  are  not  updated. 

Subroutine  OTPOT  writes  the  six  predictor  variables  and 
the  two  response  variable  values  to  a  file.  This  file  can 
then  be  read  directly  by  a  statistical  package,  such  as  SAS. 
Verification 

The  network  model  and  FORTRAN  subroutines  were  verified 
as  the  model  was  being  built.  The  SLAM  II  MONTR/TRACE 
option  allows  entities  and  attributes  to  be  monitored  all 
the  way  through  each  part  of  the  network  to  ensure  the  model 
was  operating  as  desired.  Likewise,  all  input  parameters 
were  written  to  files  to  ensure  that  the  output  made  sense 
for  that  input. 

Validation 

Because  the  real-world  system  being  modeled  has  not  yet 
reached  anything  close  to  steady  state,  it  is  difficult  to 
compare  the  model  output  to  actual  data.  And  since  the  Pave 
Low  system  is  so  small  and  unique,  it  is  not  easy  to  make 
comparisons  to  another  weapon  system  or  to  the  helicopter 
system  as  a  whole.  However,  as  the  Pave  Low  system 
approaches  steady  state,  data  will  become  available  to  help 
modify  and  validate  the  model. 
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Bxparimantation 

Bxparimantal  Design.  Examples  of  experimental  designs 
are  located  in  appendix  E.  The  first  portion  of  the 
analysis  will  consider  a  number  of  factors/variables  that 
are  suspected  as  being  significant  and  screen  out  those 
which  are  Insignificant.  With  a  large  pool  of  possibly 
significant  factors,  the  first  step  Is  to  group  homogeneous 
factors  and  use  regression  analysis  to  screen  out  entire 
groups  of  Insignificant  factors.  The  remaining  factors  are 
ungrouped  and  Individually  screened  for  significance.  The 
factors  remaining  after  this  two~stage  screening  process  are 
then  used  to  build  the  first-order  metamodel .  A  metamodel 
Is  a  model  of  a  model.  In  this  research,  It  Is  a  polynomial 
equation  that  "captures'*  what  the  simulation  model  Is  doing. 
A  lack  of  fit  test  reveals  whether  or  not  the  response 
surface  Is  a  plane  or  Is  curved,  the  latter  requiring  a 
second-order  metamodel.  This  process  of  building  a  response 
surface  Is  pictured  In  figure  4. 

In  this  particular  research,  group  screening  was  not 
needed  because  only  six  factors  were  examined.  The  factor 
screening  stage  of  the  analysis  screened  the  six  predictor 
variables  or  factors  to  determine  which  have  significant 
Impacts  on  the  two  response  variables. 
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Figure  4.  Metamodel  Building  Process 
(Bauer,  1988) 


The  six  predictor  variables  are-' 

1)  INP  (Input)  -  the  average  number  of  pilots  who  enter  the 
Pave  Low  system  each  year,  either  directly  from  flight 
school  (Undergraduate  Helicopter  Training)  or  by  cross¬ 
training  from  another  type  of  helicopter. 

2)  OHT  -  the  percentage  of  "INP”  who  come  directly  from 
flight  school,  with  no  previous  operational  flying 
experience . 

3)  TRL  (tour  length)  -  the  length,  In  months,  of  the  average 
flying  tour/asslgnment . 

4)  BLK  (block)  -  the  number  of  blocked  cockpits.  Blocked 
cockpits  are  flying  positions  reserved  for  previously 
qualified  Pave  Low  pilots  who  are  currently  serving  In  the 
rated  supplement,  but  who  eventually  return  to  the  cockpit. 
The  more  blocked  cockpits,  the  fewer  new  Inputs  can  be 
brought  Into  the  system. 

5)  PAY  (Payback)  -  the  amount  of  time  (In  flying  hours)  a 
trained  Pave  Low  pilot  Is  expected  to  remain  In  the  cockpit. 
The  greater  the  payback  requirement,  the  more  experienced 
the  pilot  force  will  be. 

6)  HRS  (hours)  -  the  average  number  of  flying  hours  each 
pilot  files  each  month.  The  more  hours,  the  more 
experienced  the  force. 

The  two  response  variables  are: 

1)  TOT  -  Total  operational  pilot  force.  This  Includes  only 
line  pilots  In  full-time,  operational  flying  (excludes 
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staff,  school  instructors,  etc.), 

2)  EXP  -  The  percentage  of  "TOT"  who  are  "experienced", 
meaning  they  have  850  total  flying  hours,  300  of  which  are 
in  the  Pave  Low  helicopter.  The  desired  operating  conditions 
for  the  model  are  defined  in  terms  of  TOT  and  EXP.  The 
number  of  authorized  line  pilots  is  approximately  100,  so 
TOT  values  that  vary  significantly  from  100  suggest  an  over 
or  under-manned  organization.  Likewise,  the  Pave  Low 
experience  level  (EXP)  should  remain  at  or  above  60%  (AFR 
51-2).  The  two  response  variables  must  be  examined 
together.  In  other  words,  it  is  not  acceptable  to  have  an 
EXP  of  70%  if  TOT  is  very  low  or  very  high. 

A  full  factorial,  two  level  experimental  design  (2  )  with 
replications  was  used.  The  Statistical  Analysis  System 
(SAS)  will  be  used  to  create  an  Analysis  of  Variance  (ANOVA) 
table  for  each  response  variable.  The  F  statistic  will  then 
be  used  to  determine  model  adequacy  and  significance  of 
factor  effects .  Furthermore ,  assumptions  regarding  constant 
variance  and  normality  of  residuals  will  be  examined  using 
residual  plots  (Montgomery,  1984:85-93,  223-224). 

From  the  regression  analysis,  a  first-order  roetamodel 
will  be  developed.  That  model  will  be  checked  for  lack  of 
fit  to  determine  if  a  second-order  model  is  required.  If 
so,  an  appropriate  three- level  design  (depending  on  the 
number  of  significant  linear  factors)  will  be  used  to 
construct  a  second-order  roetamodel. 
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Goal  Programming  Model .  Once  reasonable  response  surface 
equations  for  the  two  response  variables  have  been 
developed,  those  two  equations  will  be  set  equal  to  constant 
values  that  reflect  the  goals  of  100  pilots  with  a  60% 
experience  level.  These  goals,  with  associated  deviation 
variables,  will  become  constraints  in  a  mathematical 
programming  problem.  The  objective  function  will  seek  to 
minimize  deviations  from  the  goals.  Upper  and  lower  bounds 
will  be  placed  on  the  predictor  variables  to  restrict  them 
to  the  experimental  region .  The  solution  to  the  goal 
programming  problem  should  provide  optimal  settings  for  the 
predictor  variables  in  order  to  achieve,  as  closely  as 
possible,  the  stated  goals. 

This  chapter  described  the  Pave  Low  model  and  the 
experimental  design  for  the  output  analysis.  The  next 
chapter  discusses  the  results  of  that  output  analysis . 
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IV.  Analyaia 


Introduction 

The  purpose  of  this  chapter  Is  to  describe  the  operation 
of  the  SLAM  model  and  the  analysis  of  the  model  output  using 
response  surface  methodology  and  goal  programming.  The 
response  surface  methodology  Is  used  to  build  response 
equations  that  are  then  used  as  goals  In  the  goal 
programming  formulation.  The  goal  programming  model,  In 
turn,  determines  optimal  values  for  the  decision  variables. 
SLAM  II  Model  Operation 

To  estimate  the  length  of  the  transient  period,  several 
model  runs  were  conducted  for  time  periods  of  100  years. 

The  time  average  of  the  experience  level  was  then  plotted 
(See  appendix  G) .  The  plot  indicates  the  model  has  a  warm¬ 
up  period  of  approximately  150  months  before  reaching  a 
steady  state  condition.  In  subsequent  runs,  all 
statistical  arrays  were  cleared  after  200  months .  The  model 
was  then  allowed  to  run  for  another  700  months  or  about  60 
years . 

Output  Analyaia 

Factor  Screening.  The  model  was  run  using  ten 
replications  of  the  six-factor,  two-level  full  factorial 
design  with  the  six  predictor  variables  shown  in  table  2. 
These  variables,  along  with  the  two  response  variables  were 
analyzed  using  SAS  procedures  PROG  REG  and  PROG  STEPWISE. 
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Table  2 .  Factors  and  Settings 


EACTQB 

LOW 

HIGH 

I  NT 

12 

24 

OHT 

0 

.6 

TRL 

36 

60 

PAY 

500 

750 

6LK 

0 

30 

HRS 

10 

15 

The  residual  plot  from  the  second  response  variable,  TOT, 
Indicated  a  slightly  increasing  variance.  A  natural  log 
transformation  of  the  response  corrected  the  problem  (All 
residual  plots  and  normality  plots  are  in  appendix  F) .  The 
resulting  ANOVA  tables  are  in  tables  3  and  4. 


^  Table  3.  ANOVA  for  EXP  factor  screening 


HR 

SOM  OF 
SOOARES 

MEAN 

SQOARE 

HHRI 

MODEL 

6 

93673.59 

15612.26 

866.10 

0.0001 

INP 

1 

758.21 

758.21 

42.06 

0.0001 

OHT 

1 

65642.40 

65642.40 

3641.54 

0.0001 

TRL 

1 

3691.20 

3691.20 

204.77 

0.0001 

BLK  . 

1 

13987.60 

13987.60 

775.97 

0.0001 

PAY 

1 

6347.88 

6347.88 

352.15 

0.0001 

HRS 

1 

3246.30 

3246.30 

180.09 

0.0001 

ERROR 

633 

11410.45 

18.026 

TOTAL 

639 

105084.04 

ADJ 

R"^  =  .89 
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Table  4.  ANOVA  for  TOT  factor  screening 


SOM  OF 
SQUARES 

MEAN 

SQUARE 

RIPIVRili 

MODEL 

6 

98.38 

16.40 

1006.49 

0.0001 

INP 

1 

60.89 

60.89 

3737.30 

0.0001 

OHT 

1 

6.64 

6.64 

407.53 

0.0001 

TRL 

1 

2.68 

2.68 

164.32 

0.0001 

6LK 

1 

11.34 

11.34 

695.92 

0.0001 

PAY 

1 

5.92 

5.92 

363.17 

0.0001 

HRS 

1 

10.93 

10.93 

670.72 

0.0001 

ERROR 

TOTAL 

633 

639 

10.31 

108.70 

.016 

ADJ 

R^  =  .90 

The  p  values 

of  0.0001 

for  the  six 

predictor 

variables 

indicate  that  all  six  are  significant.  It  nas  noted, 
however,  that  UHT  was,  by  far,  the  most  significant  factor 
to  EXP,  and  INP  was  the  most  significant  to  TOT  (based  on 
significantly  larger  F  values).  To  further  simplify  the 
linear  metamodel,  then,  the  remaining  four  factors  were 
screened  out. 

First  Order  Metamodal  Ten  replications  of  a  three- 
level,  two-factor  full  factorial  experimental  design  were 
run.  A  two-level,  two-factor  design  was  used  in  a  trial  run 
and  Indicated  that  a  second  order  model  was  needed.  The 
four  screened  factors  were  set  to  constant,  nominal  values 
so  they  would  not  Influence  the  results.  Factor  settings 
are  shown  in  table  5 . 

The  two-factor  linear  model  residual  plots  were  convex 
curves,  indicating  quadratic  tendencies,  which  were 
confirmed  by  using  SAS  procedure  PROC  RSRE6.  This  procedure 
breaks  the  regression  or  model  sum  of  squares  into  linear. 


35 


Table  5.  Factor  Settings  for  First-Order  tietamodel 


FAOTOR 

LOW 

ttUL 

HIGH 

IHP 

12 

16 

24 

DHT 

0 

.25 

.50 

TRL 

48 

PAY 

680 

BLK 

15 

HRS 

12.5 

quadratic,  and  cross-product  components.  Tables  6  and  7 
show  both  significant  quadratic  and  cross-product 
components.  If  a  two-level  design  had  been  used  instead  of 
a  three-level  design,  PROG  RSREG's  lack  of  fit  test  would 
have  indicated  a  lack-of-fit  for  a  linear  model .  That  was 
the  case  in  the  trial  run. 


Table  6.  First-order  lack-of-fit  for  EXP 


m^RURiii 

SUM  OF 
SQUARES 

MEAN 

SQUARE 

LINEAR 

2 

13256.03 

6628.41 

18725.26 

0.0001 

QUADR 

2 

25.69 

12.84 

36.29 

0.0001 

CROSS 

1 

7.26 

7.26 

20.51 

0.0001 

TOTAL 

5 

13289.78 

2657.96 

7508.72 

0.0001 

RESIDUAL 

L.O.F. 

3 

2.65 

0.88 

2.57 

0.0552 

PURE  E. 

171 

58.94 

0.34 

o 

TOTAL 

174 

61.59 

0.35 

ADJ  R 

=  .99 

Table  7 .  First-order  lack-of-f it  for  TOl 


SUM  OF 
SQUARES 

MEAN 

SQUARE 

LINEAR 

2 

165973.55 

82986.78 

62480.93 

0.0001 

QOADR 

2 

5739.51 

2869.76 

2160.64 

0.0001 

CROSS 

1 

753.38 

753.38 

567.22 

0.0001 

TOTAL 

5 

172466.44 

34493.29 

25970.07 

0.0001 

RESIDUAL 

L.O.F. 

3 

18.70 

6.23 

5.02 

0.0025 

PURE  E. 

171 

212.41 

1.41 

0 

TOTAL 

174 

231.11 

1.33 

ADJ  R 

=  .99 

Second  Order  Metamodel.  The  residual  plots  for  the 
second-  order  model  based  on  the  two  factors,  INP  and  UHT, 
Indicate  a  non-constant,  decreasing  variance  for  response 
variable  EXP  and  a  non-constant,  wave  pattern  for  TOT. 

These  could  indicate  the  need  for  a  transformation.  Several 
transformations  were  tried,  including  natural  log,  arcs in, 
and  square  root  functions  of  the  response  variable .  Based 
on  residual  plots,  the  best  transformations  were  ln(84-EXP) 
and  In (TOT).  ANOVA  tables  are  in  tables  8  and  9. 


Table  8.  ANOVA  for  second-order  metamodel  (EXP) 


RH 

SUM  OF 
SQUARES 

MEAN 

SQUARE 

MODEL 

5 

71.48 

14.30 

8770.87 

0.0001 

INP 

1 

1.27 

1.27 

784.38 

0.0001 

UHT 

1 

66.33 

66.33 

40913.57 

0.0001 

DHT*UHT 

1 

3.63 

3.63 

2241.22 

0.0001 

INP*UHT 

1 

0.25 

0.25 

152.70 

0.0001 

ERROR 

TOTAL 

175 

179 

0.28 

71.77 

0.0016 

ADJ 

R^  =  .99 
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Table  9.  ANOVA  for  second-order  metamodel  (TOT) 


SUM  OF 
SQUARES 

MEAN 

SQUARE 

l^fliRllliiil 

MODEL 

5 

11.63 

2.33 

25580.85 

0.0001 

INP 

1 

10.87 

10.87 

120045.03 

0.0001 

DHT 

1 

0.63 

0.63 

6968.06 

0.0001 

INP* INP 

1 

0.13 

0.13 

1450.75 

0.0001 

INP*0HT 

1 

0.01 

0.01 

13.49 

0.0003 

ERROR 

175 

0.02 

0.00009 

TOTAL 

179 

11.65 

ADJ  R 

^  =  .99 

After  decoding,  the  second  order  metamodel  for  response 
EXP  was: 


ln(84-EXP)  =  1.3312  +  0.0265*INP  +  6.0524*0HT 
(.007)  (.004)  (.004) 

-  4.8224*0HT^  -  0 . 0371*INP*UHT 
(.006)  (.006) 


The  second-order  metamodel  for  TOT  was: 

In(TOT)  =  4.2696  -  0.0078»INP  +  0.2432*UHT 

(.002)  (.001)  (.001) 

+  0.0016*INP^  +  0.0026*INP*DHT 

(.002)  (.001) 


Figure  4  shows  both  responses  plotted  In  the  INP/UHT 
plane.  Since  both  response  surfaces  were  second  order,  they 


are  both  paraboloids.  The  figure  is  a  plot  of  the 


projection  of  the  paraboloid  contours  in  the  INP/OHT  plane. 
The  contours  are  nearly  straight  lines  because  the  region  of 


Interest  is  so  far  from  the  parabolas'  vertices.  This 


figure  will  vary  with  the  values  used  for  the  other  four 
decision  variables  (TRL,  PAY,  BLK,  AND  HRS).  This 
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Goal  programming  Is  a  form  of 


mathematical  programming  that  allows  more  than  one  objective 
to  be  achieved  In  the  same  problem.  So  rather  than 
minimizing  or  maximizing  one  objective,  goal  programming 
minimizes  deviations  from  one  or  more  established  goals. 
Goals  C2ua  be  prioritized  and  deviations  can  be  weighted 
depending  on  the  judgment  of  a  decision  maker  or  expert 
(Hllller  and  Lleberman,  1986:  242-246). 

The  two  goals  In  this  research  are:  1)  to  maintain  the 
force  (TOT)  at  close  to  100  pilots  and  2)  to  maintain  the 
force  experience  level  (EXP)  at  a  minimum  of  60%.  A  linear 
goal  programming  problem  was  formulated  as  follows,  using 
the  six-factor  linear  metamodels  In  the  goal  equations: 
Minimize:  dj^~  +  dg”  +  dg^ 

Subject  to:  1)  36.79  -  0.81INP  -  40.52DHT  +  0.2TRL  + 

0.03PAY  +  .31BLK  +  0.90HRS  +  d^~  -  dj^"^  =  60 

2)  -19.78  +  6.05INP  +  45.36UHT  +  0.64TRL  + 

0.09PAY  +  0.92BLK  -  6.20HRS  +  dg"  -  dg'*'  =  100 

3)  12  <  INP  <  24 

4)  0  <  OHT  <  .5 

5)  36  <  TRL  <  60 

6)  500  <  PAY  <  750 

7)  0  <  BLK  <  30 

8)  10  <  HRS  <  15 

The  "d"  variables  are  the  deviation  variables  for  the  two 
goals.  Since  the  experience  level  can  be  allowed  to  go 
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above,  but  not  below  60X,  was  left  out  of  the  objective 

function.  The  remaining  three  deviation  variables  are 
minlffllsed.  All  the  predictor  variables  have  lower  and  upper 
bounds,  to  keep  the  solution  set  within  the  experimental 
region . 

The  goal  programming  model  was  run  on  LIHDO  (Schrage, 
1986),  a  linear  programming  package.  The  results,  after  13 
iterations ,  were : 

INP  =  17.92 
OHT  =  0.12 

TRL  =  36.0 

PAY  =  500.0 
BLK  =  0.0 

HRS  =  10.0 

Objective  function  value  =  0.0 
All  deviation  variables  =  0,0 

Thus,  the  goals  were  able  to  be  met  exactly  by  setting 
all  but  the  first  two  predictor  variables  to  their  lower 
bound.  These  optimal  values  were  run  in  the  Pave  Low 
simulation  model  and  produced  an  EXP  value  of  approximately 
61  percent  and  a  TOT  value  of  109.  While  not  exactly  on 
target,  they  are  an  Indication  the  goal  programming  model 
accurately  reflects  the  simulation  model  output. 

In  contrast,  the  second-order  metamodels  were  used  as 
goals  in  a  non-linear  goal  programming  problem.  The  optimal 
solution  was  to  set  INP  to  14.2  and  DHT  to  42  percent.  When 
this  solution  was  run  in  the  simulation  model,  it  met  the 
experience  goal,  but  exceeded  the  manning  (TOT)  goal  by 
about  22.  The  difference  beteen  the  two  metamodels  (linear 
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vs  quadratic)  is  the  linear  includes  the  four  less 
significant  decision  variables.  These  four  variables  help 
provide  a  better  reflection  of  the  simulation  in  the 
metamodel  equations. 

The  difficulty  with  any  "optimal”  solution  is  that  it 
may  not  be  practical.  Although  all  the  predictor  variables 
are  decision  variables  and  can  be  controlled,  some  are  more 
easily  varied  than  others.  Since  INP  and  OHT  are  probably 
the  easiest  to  change,  the  other  predictor  variables  can  be 
set  to  values  that  seem  to  most  closely  reflect  reality  and 
let  the  model  solve  for  IHP  and  UHT. 

This  chapter  discussed  the  analysis  of  the  simulation 
output.  In  light  of  this  analysis,  the  next  chapter  will 
make  policy  recommendations  as  well  as  recommendations  for 
further  research. 
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This  chapter  summarizes  the  research  and  its 
Implications,  and  then  provides  recommendations  for  policy 
decisions  and  further  research. 

Summary  of  Reaearch 

The  objective  of  this  research  was  to  model  the  Pave  Low 
pilot  career  field  and  to  examine  how  best  to  sustain  an 
experienced  pilot  force.  This  was  accomplished  by 
developing  a  simulation  model  that  permitted  experimentation 
with  key  policy/decision  variables.  Response  surface 
methodology  was  then  used  to  screen  out  less  Important 
variables  and  to  develop  response  equations  which  captured 
the  simulation  outcomes  in  simple  polynomials.  These 
polynomials  were  used  to  define  goals  within  a  goal 
programming  model.  That  model  produced  optimal  values  for 
the  policy  variables  in  order  to  achieve  the  desired  manning 
and  experience  levels  for  the  Pave  Low  pilot  force. 
Implicationa 

The  value  of  this  research  is  in  developing  a  methodology 
both  to  describe  the  Pave  Low  pilot  career  field  and  to 
determine  optimal  policy  decisions.  As  more  complete  data 
on  the  pilot  force  becomes  available,  the  models  used  in 
this  research  can  be  easily  modified  to  more  accurately 
reflect  reality. 
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This  research  suggests  two  policy  recommendations  once 


the  Pave  Low  pilot  force  Is  fully  manned.  First,  only 
about  ten  percent  of  the  pilots  entering  the  Pave  Low  should 
be  right  out  of  Undergraduate  Helicopter  Training  (UHT) 
to  avoid  causing  the  experience  level  to  drop  below  60 
percent.  Second,  the  Input  rate  Into  the  system  will 
probably  have  to  be  reduced  slightly  from  the  current  twenty 
pilots  per  year  to  prevent  over-manning.  The  exact  amount 
of  reduction  will  depend  on  future  retention  rates . 
Recommendation  for  Future  Research 

Future  research  should  focus  on  the  Implications  of 
overseas  basing,  such  as  how  short  overseas  tours  and  the 
ratio  of  overseas  to  stateside  assignments  affects  manning 
and  experience  levels . 


Appendix  A :  Rfltflntion  Rates 


Year  of  Service 

1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 


Retention  Rate 

1.000 

.986 

1.000 

.958 

.974 

.971 

.909 

.963 

.955 

.977 

.875 

.886 

.892 

1.000 

1.000 

.983 

1.000 

.955 

.906 

.745 

.521 

.667 

1.000 

.778 

.500 


Appendix  B: 


SLAB  11  -HatHork  Diagram 


Appendix  C: 


GEN, MCKOY.PAVEl, 6/1 5/1988, 32, N.N. ,N,N,72; 

LIMITS, 5, 8, 750; 

TIMST, XX( 4) , PERCENT  EXPERIENCED; ,20/0./. 05; 
TIMST,XX(1) ,NO  OF  EXP  PILOTS; 

TIMST, XX(2) , NO  OF  ACTIVE  PILOTS; 

SEEDS, 1145661099(1); 

; RECORD, TNOW, TIME, ,P,5.0; 

;VAR,XX(4), . , PERCENT  EXPERIEN, 0. 0, 1 . 0; 
EQDIVALENCE/XX( 1 ) , EXPD/ 

XX(2),ACTV/ 

XX(3),REQD/ 

XX(4),PERC/ 

XX{5) ,EXCS; 

NETWORK; 

; CLEAR  STATISTICAL  ARRAYS  AT  TNOW=200 
CREATE, ,200; 

CLEAR  EVENT, 3; 

TERM; 

; CREATE  A  CLASS  OF  2  STUDENTS  W/  T.B.C.  OF  XX(6) 
CREATE,XX(6); 

ASSIGN,ATRIB(8)=2,1; 

0NBATCH,8,1; 

QDEl  QDEOE(l); 

ACT(20); 

; LABEL  STUDENTS  AS  UHTS  OR  PQPS 
GOON,l; 

ACT/1, , XX (7), UHT; UHTS 
ACT/2 ,, 1-XX( 7 ), POP ; PQPS 

;UHT  IS  ASSIGNED  TIME  IN  SERVICE  OF  12  MONTHS 
UHT  ASSIGN , ATRIB( 3 ) =0 . , 

ATRIB(4)=12. , 

ATRIB(5)=0., 

ATRIB(7)=0. ,1; 

ACT, , ,G1; 

;PQP  IS  ASSIGNED  T.I.S.  lAW  TRIANGULAR  DISTRIBUTION 
PQP  ASSIGN, ATRIB( 3 ) =1. , 

ATRIB(4)=TRIAG(42. ,48. ,54. ,1), 
ATRIB(5)=0. , 

ATRIB(7)=0. ,1; 

;CCTS  TAKES  7  MONTHS, 80  FLYING  HOURS 
G1  GOON,l; 

ACT/3 , 7 . ; CCTS 
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ASSIGN,ATRIB(4)=ATRIB(4)-»-7. , 

ATRIB(5)=ATRIB(5)+80. , 

ACTV=ACTV+1 . , 

PERC=EXPD/ACTV , 1 ; 

; STUDENT  IS  NOW  AN  ACTIVE,  THOUGH  INEXPERIENCED  PILOT 

; PILOT  GETS  FIRST  OPERATIONAL  ASSIGNMENT 
;TOOR  LENGTH  I AW  TRIANGULAR  DISTRIBUTION 
QUE2  QUEUE(2); 

ACT(IOO) ; 

ASSIGN, ATRIB(1)=TRIAG(XX(8),XX(9),XX(10),1),1; 

;SEE  FORTRAN  CODE 
El  EVENT, 1,1; 

;OHT  GETS  A  2ND  TOUR  TO  GET  850  HRS 

ACT/4, ,ATRIB(3) .EQ.0.0.AND.ATRIB{5) .LT.850. ,G2;GET 

EXPER 

;DHTS  W/  850  HRS 

ACT/6 , , ATRIB( 3 ) . EQ . 0 . , UHTX ; EXP-UHT 
;PQPS  WHO  ATTRIT  FROM  SYSTEM  W/0  GETTING  300  HRS 

ACT/8, ,ATRIB(5) .LE. 300. , LOSE, -FIRST  YR  LOSS/PQP 
;PQPS  W/  300  HRS 

ACT/10, ,ATRIB(3) .EQ.l. ,PQPX; EXP-POP 

G2  GOON,l; 

ACT/12, ATRIB( 1 ), ,G3; 

G3  GOON.l; 

;0HTS  WHO  ATTRIT  DURING  FIRST  TOUR 

ACT/14, ,ATRIB(7) . EQ. 1 . 0 , LOSE; SEPARl 
ACT/16, , ,QUE2; 

; COMPUTES  #OF  MONTHS  INTO  TOUR  THAT  PILOT 
; REACHES  EXPERIENCED  STATUS  (850  HRS) 

UHTX  ASSIGN, EXCS=ATRIB( 5) -850. , 

REQD=ATRIB(2)-EXCS, 

ATRIB(6)=REQD/15. ,1; 

;0HTS  WHO  ATTRIT  DURING  2ND  TOUR  W/0  REACHING  850  HRS 
ACT/18, ,ATRIB( 6) .GE.ATRIB(l) ,L0; 

;DHTS  W/  850  HRS 

ACT/20, ATRIB(6) ,ATRIB( 6) .LT.ATRIB(l) ,CNTR;#EXP  DHTS 

LO  GOON , 1 ; 

ACT/22, ATRIB(l) , ,G4; 

; COMPOTES  *  OF  MONTHS  INTO  TOUR  THAT  PILOT 
; REACHES  EXPERIENCED  STATUS  (300  HRS) 

PQPX  ASSIGN , EXCS=ATRIB( 5 ) -300 . , 

REQD=ATRIB(2)-EXCS, 

ATRIB(6)=REQD/15. ,1; 
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;PQPS  WHO  ATTRIT 

ACT/24, ,ATRIB(6) .GE.ATRIB(1),L1;PQP  NEVER  EXP 
;PQPS  W/  300  HRS 

ACT/26 , ATRIB ( 6 ) . ATRIB ( 6 ) . LT . ATRIB ( 1 ) , CNTR ; #EXP  PQPS 

; COUNTS  EXPERIENCED  PILOTS,  COMPOTES  PERCENT  EXP 
CNTR  ASSIGN ,  EXPD=EXPD->-l . , 

PERC=EXPD/ACTV, 

; COMPUTES  TIME  REMAINING  IN  TOUR 

ATRIB( 1 ) =ATRIB{ 1 ) -ATRIB( 6 ) , 

ATRIB(6)=0.0,1; 

ACT/28, ATRIB( 1) , ,G4; TOTAL  EXP 

LI  GOON.l; 

ACT/3 0, ATRIBd ), ,G4;PQP  ONLY 
G4  GOON.l; 

.•EXPERIENCED  PILOTS  WHO  ATTRIT 

ACT/32, ,ATRIB(7) .EQ. 1.0. AND. ATRIB{6) .EQ. 0.0, DEPT; 

; NEVER  EXPERIENCED  PILOTS  WHO  ATTRIT 

ACT/34, ,ATRIB(7) .EQ. 1.0. AND. ATRIB(6) .NE. 0.0, LOSE; 

; EXPERIENCED  UHTS  WHO  OWE  PAYBACK  FOR  TRAINING 

ACT/36, ,ATRIB(3) .EQ.0.0.AND.ATRIB(5) .LT. 1000. ,PAY0; 

: EXPERIENCED  PQPS  WHO  OWE  PAYBACK  FOR  TRAINING 

ACT/38, ,ATRIB(3).EQ. 1.0. AND. ATRIB(6) .LT.XX( 11), PAYl; 
; PILOTS  W/  MORE  THAN  11  YRS  SERVICE 

ACT/40, ,ATRIB(4) .GE . 132 . ,QDE4; SKIP  JOB3 
; EVERYONE  ELSE 

ACT/42. , ,Q0E3; TO  JOB2 

; ASSIGN  SUBSEQUENT  FLYING  TOUR  LENGTHS  TO  UHTS 
PAYO  ASSIGN, ATRIB(1)=TRIAG(XX(8), XX(9) ,XX(10) ,1) ,1; 

EVENT ,1,1; 

ACT , , , LO ; 

.•ASSIGN  SUBSEQUENT  FLYING  TOUR  LENGTHS  TO  PQPS 
PAYl  ASSIGN,ATRIB(1)=TRIAG(XX(8),XX(9) .XX(IO) ,1) ,1; 

EVENT, 1,1; 

ACT, , ,L1; 

;7  TO  11  YEAR  GROUP 
QUE3  QUEUE(3); 

ACT(100)/44;TO  G5 
G6  GOON.l; 

;n  OF  BLOCKED  COCKPITS  FOR  7-11  YR  GROUP 

ACT, ,NNACT(46) .LT.XX(12) ,FLY2;FLY2 
; OTHERWISE  TO  NON-FYING  JOB 
ACT , , , NOF2 ; NOFLY2 


; FLYING  ASSIGNMENTS  FOR  7-11  YR  GROUP 

FLY2  ASSIGN, ATRIB( 1 )=TRIAG(XX{8), XX(9) .XX( 10) .1),1; 
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EVENT, 1,1; 

ACT/46 , ATRIB( 1) , , G6 ; FLY2 

; NON-FLYING  ASSIGNMENTS  FOR  7-11  YR  GROUP 

NOF2  ASSIGN , ATRIB( 1 ) =TRIAG( 30 . , 36 . , 42 . , 1 ) , 

ACTV=ACTV-1 . , 

EXPD=EXPD-1 . , 

PERC=EXPD/ACTV , 1 ; 

EVENT ,2,1; 

ACT/48 , ATRIB( 1 ) ; NOFLY2 
ASSIGN , ACTV=ACTV+1 . , 

EXPD=EXPD+1 . , 

PERC=EXPD/ACTV,1; 

G6  GOON,l; 

; ATTRITION 

ACT/50, ,ATRIB(7) . EQ . 1 . 0 , DEPT ; SEP AR3 

;LESS  THAN  11  YRS  TIME  IN  SERVICE 

ACT/52, ,ATRIB(4) .LT. 132. ,Q0E3;<11YRS  GATE  TIME 

; EVERYONE  ELSE 

ACT/54; TO  JOB  3 

;  12-17  YR  GROUP 

QUE4  QUEUE(4); 

ACT(120)/56;TO  G7 

G7  GOON,l; 

; BLOCKED  COCKPITS  FOR  12-17  YR  GROUP 

ACT, ,NNACT(58) .LT.XX{12),FLY3;FLY3 

ACT, , ,NOF3;NOFLY3 

FLY3  ASSIGN , ATRIB( 1 ) =TRIAG ( XX( 8 ) , XX( 9 ) , XX( 10 ) , 1 ) , 1 ; 
EVENT, 1,1; 

ACT/58 , ATRIB( 1 ) , , G8 ; FLY3 

NOF3  ASSIGN,ATRIB(1)=TRIAG(30. ,36. ,42. ,1), 
ACTV=ACTV-1. , 

EXPD=EXPD-1 . , 

PERC=EXPD/ACTV,1; 

EVENT, 2,1; 

ACT/60 , ATRIB( 1 ) ; NOFLY3 
ASSIGN , ACTV=ACTV-H  . , 

EXPD=EXPD+1 . , 

PERC=EXPD/ACTV,1; 

G8  GOON , 1 ; 

ACT/62 , , ATRIB( 7 ) . EQ . 1 . 0 , DEPT ; SEPAR4 

ACT/64, ,ATRIB(4) .LT.204. ,QDE4;<17YRS 

ACT/66 ; LAST JOB 


;  18-20+  YR  GROUP 
QDE5  Q0EUE(5); 

ACT(160)/68;TO  G9 
G9  G00H,1; 

BLOCKED  COCKPITS  FOR  SENIOR  OFFICERS 
ACT , , NNACT( 70 ) . LT . 0 , FLY4 ; FLY4 
ACT, , ,N0F4;N0FLY4 

FLY4  ASSIGN,ATRIB(1)=TRIAG(18. ,24. ,30. ,1) ,1; 

EVENT, 1,1; 

ACT/70,ATRIB(1) , ,G10;FLY4 

NOF4  ASSIGN,ATRIB(1)=TRIAG(30. ,36. ,42. ,1), 

ACTV=ACTV-1 . , 

KXPD=EXPD-1 . , 

PERC=EXPD/ACTV,1; 

EVENT, 2,1; 

ACT/72 , ATRIB( 1 ) ;N0FLY4 
ASSIGN, ACTV=ACTV+1 . , 

EXPD=EXPD+1 . , 

PERC=EXPD/ACTV,1; 

GIO  GOON.l; 

ACT/74, ,ATRIB(7) .EQ.1.0.AND.ATRIB(4) .LT.240. ,DEPT; 
ACT/76, ,ATRIB(4) .LT.240. ,QDE5; <20YRS 
; RETIREES 

ACT/78; RETIRE 
ASSIGN,EXPD=EXPD-1. , 

ACTV=ACTV-1 . , 

PERC=EXPD/ACTV,1; 

COLCT , ALL , RETIREES , . 1 ; 

TERM; 

;  PILOTS  LOST  BEFORE  BECOMING  EXPERIENCED 
LOSE  ASSIGN , ACTV=ACTV-1 . , 

PERC=EXPD/ACTV,1; 

ACT/80 , , , COLL ; NEVER  EXPD 

; EXPERIENCED  PILOTS  WHO  SEPARATE 
DEPT  ASSIGN , EXPD=EXPD-1 . , 

ACTV=ACTV-1. , 

PERC=EXPD/ACTV,1; 

COLL  COLCT , ALL , SEPARATED , , 1 ; 

TERM; 

END; 

INIT, 0,900; 

SIMULATE; 

;THIS  IS  JUST  ONE  LINE  OF  THE  DESIGN  MATRIX 

-1  -1  -1  -1  -1 

FIN; 
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Appendix  D :  FORTRAN  Code 


DIMENSION  NSET( 20000) 

C0MM0N/SC0M1/ATRIB( 100 ) , DD( 100 ) , DDL( 1 00 ) , DTNOW , II , MFA , 
IMSTOP , NCLNR , NCRDR , NPRNT , NNRON . NNSET , NT APE , SS ( 1 00 ) , 
2SSL(100),TNEXT,TNOW.XX(100) 

COMMON  QSET(20000) 

COMMON/DCOMl/RET ( 25 ) 

EQUIVALENCE ( NSET ( 1 ) >  QSET ( 1 ) ) 

NNSET=20000 

NCRDR=5 

NPRNT=6 

NTAPE=7 

OPEN( 1 1 , STATUS= ' NEW ' , FILE= ' THESIS . DAT ' 

,FORM=' FORMATTED') 

OPEN ( 12 , STATDS= ' NEW ' , FILE= ' SAS . DAT ' 

,FORM=' FORMATTED') 

CALL  SLAM 

STOP 

END 


SUBROUTINE  INTLC 

COMMON/SCOMl/ATRIB( 1 00 ) , DD( 100 ) , DDL ( 1 00 ) , DTNOW , I I , MFA , 
IMSTOP , NCLNR , NCRDR , NPRNT , NNRON , NNSET , NT APE , SS ( 100 ) . 
2SSL ( 1 0 0 ) , TNEXT , TNOW , XX ( 1 0  0 ) 

COMMON  QSET (20000) 

COMMON/UCOMl/RET ( 25 ) , XINP , UHT , TRL , PAY , BLK , HRS 
EQUIVALENCE  ( XX(1 ) , EXPD) , ( XX( 2 ) , ACTV ) 

READ  AND  INITIALIZE  THE  VARIABLES 
READ(NCRDR,*)XX(6) ,XX(7),XX(9) ,XX(11) .XX(12) 

C  CODED  VARIABLES 

XINP=XX(6) 

DHT=XX(7) 

TRL=XX(9) 

PAY=XX(11) 

BLK=XX(12) 

HRS=-1. 

C  ONCODED  FOR  SLAM 

XX(6)=-.5*XX(6)+1.5 

XX(7)=.25*XX{7)+.25 

XX(9)=12.*XX(9)+48. 

XX(8)=XX{9)-6.0 

XX(10)=XX(9)+6.0 

XX(ll)=125.*XX(ll)+625. 

XX(12)=(15.*XX(12)+15. )/2. 

XX(13)=10. 
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C  SET  OP  RETENTION  RATES  FOR  YEAR  OF  SERVICE  1-20 
DATA 

RET/1.0, .986,1.0, .958, .974, .971, . 909 , . 963 , . 955 

1. . 977. .875. .886. .892.1.0.1.0. .983.1.0. .955. .906. .745. 

2.521. . 667.1.0. .778. .5/ 

RETURN 

END 


SOBROOTINE  EVENT (I) 

COMMON/SCOM1/ATRIB(100) ,DD(100) ,DDL(100) ,DTNOW, II 
1 , MFA , HSTOP , NCLNR , NCRDR , NPRNT , NNRON . NNSET , NT APE 
2,SS(100) ,SSL(100) ,TNEXT,TNOW,XX(100) 

COMMON  QSET( 20000) 

COMMON/OCOMl/RET( 25 ) , NBl , NB2 , NB3 

GO  TO  (1,2,3)  I 

C  FOR  FLYING  ASSIGNMENTS 

1  CALL  FLYER 
RETURN 

C  FOR  NON-FLYING  ASSIGNMENTS 

2  CALL  NONFLY 
RETURN 

C  TO  CLEAR  STAT  ARRAYS  AFTER  200  MONTHS 

3  CALL  CLEAR 
RETURN 
END 


SUBROUTINE  FLIER 

COMMON/SCOMl/ATRIB { 100 ) , DD( 100 ) , DDL ( 1 00 ) , DTNOW , I I 
1 , MFA , MSTOP , NCLNR , NCRDR , NPRNT . NNRON , NNSET , NT APE , 
2SS(100) ,SSL(100) ,TNEXT,TNOW,XX(100) 

COMMON  QSET( 20000) 

COMMON/OCOMl/RET ( 25 ) , NBl , NB2 . NB3 

C  TIME  IN  SERVICE  IN  YEARS 
ITIS=INT(ATRIB(4)/12. ) 

C  ASSIGNED  TOUR  LENGTH  IN  YEARS 

ITRL=INT{ ATRIB( 1 )/12 . ) 

C  DRAW  FROM  UNIFORM  DISTRIBUTION 

PROB=ONFRM( 0.0, 1.0,1) 

XRBT=1.0 

ATRIB(7)=0.0 

C  FOR  AS  MANY  YEARS  AS  TOUR,  CHECK  RETENTION 
DO  10  I=1.ITRL 
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XRET=RET(ITIS+I) 

IF  TIME  IH  SERVICE  EXCEEDS  25  YEARS,  USE  YR  25 
RETENTION 
TEMP=ITIS+I 

IP ( TEMP . QT . 25 ) XRET=RET( 25 ) 

IF  SEPARATING,  COMPOTE  FIGURES  FOR  SHORTENED  TOUR 
LENGTH 

IF(PROB.GT.XRET)  THEN 
ATRIB(1)=I*12. 

ATRIB( 4 ) =ATRIB( 4) +ATRIB( 1 ) 

ATRIB ( 2 ) = 1 0 . *ATRIB ( 1 ) 

ATRIB ( 5 ) = ATRIB( 5 ) +ATRIB ( 2 ) 

ATRIB(7)=1.0 
GO  TO  3 
ENDIF 

C  IF  NOT  SEPARATING,  COMPUTE  FIGURES  FOR  ASSIGNED  TOUR 
10  CONTINUE 

ATRIB ( 2 ) = 1 0 . *ATRIB ( 1 ) 

ATRIB ( 4 ) = ATRIB ( 4 ) +ATRIB ( 1 ) 

ATRIB( 5 ) =ATRIB( 5 ) +ATRIB( 2 ) 

3  RETURN 
END 


SUBROUTINE  NONFLY 

COMMON/SCOMl/ATRIB( 1 00  > , DD( 1 00) . DDL( 100 ) , DTNOW , 

1 I I , MFA , MSTOP , NCLNR , NCRDR , NPRNT , NNRUN , NNSET , NT APE 
2,SS(100),SSL(100) ,TNEXT.TNOW,XX(100) 

COMMON  QSET( 20000) 

COMMON/OCOMl/RET ( 25 ) , NBl , NB2 , NB3 

ITIS=INT(ATRIB(4)/12.) 

ITRL=:INT(ATRIB(1)/12. ) 

PROB=UNFRM( 0.0, 1.0,1) 

XRET=1.0 
ATRIB(7)=0.0 
DO  10  I=1,ITRL 
XRET=RET(ITIS+I) 

TEMP=ITIS+I 

IF(TEMP.GT.25)XRET=RET(26) 

IF(PROB.GT.XRET)  THEN 
ATRIB(1)=I*12. 

ATRIB( 4 ) =ATRIB( 4 ) +ATRIB( 1 ) 

ATRIB(7)=1.0 
GO  TO  3 
ENDIF 

10  CONTINUE 

ATRIB(4)=ATRIB(4)+ATRIB{1) 

3  RETURN 
END 
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SDBRODTINE  OTPUT 

COMMON/SCOMl/ATRIBdOO)  ,DD(100)  ,DDL(  100)  .DTNOW 
1,11, UFA , MSTOP , NCLNR , NCRDR , NPRHT , NNRON , NNSET , NT APE 
2,SS(100) .SSL(IOO) .TNEXT,TNOW,XX(100) 

COMMON  QSET( 20000) 

COMMON/OCOMl/RET ( 25 ) , XINP , DHT , TRL , PAY , BLK , HRS 
EQUIVALENCE  (XX(4) ,PERC) , (XX(2) ,ACTV) 
TEMP=100*TTAVG(1) 

C  WRITE  ONCODED  VARIABLES  AND  OOTPOT  TO  A  FILE 

WRITE(11,100)XX(6),XX(7),XX(9).XX(11).XX(12).XX(13) 

1,TBMP,TTAVG(3) 

C  WRITE  CODED  VARIABLES  AND  OOTPOT  TO  A  FILE 

WRITE ( 1 2 , 20 0 ) XINP , OHT , TRL , PAY , BLK , HRS , TEMP , TT AVG ( 3 ) 

100  F0RMAT(8(1X,F6.2,1X)) 

200  FORMAT(6(2X,F5.2.1X) ,2(1X,F5.1.1X)) 

RETORN 

END 
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Factor  Screen :  EXP 
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Normal  Probability  Plot:  EXP 
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Hormal  Probability  Plot:  In (TOT) 
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First-Order  Metamodel :  TOT 
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Second-Order  Metamodel:  ln( 84-EXP) 
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Hormal  Probability  Plot:  ln( 84-EXP) 
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Second-Order  Metamodel:  In (TOT) 
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Appendix  G:  Steady-State  Graph 
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